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ABSTRACT 

Strong observational evidence for a fluctuating ultraviolet background (UVB) has 
been accumulating through a number of studies of the H i and He n Lya forest as 
well as accurate IGM metallicity measurements. UVB fluctuations could arise both 
from the inhomogeneous distribution of the ionizing sources and/or from radiative 
transfer (RT) through the filamentary IGM. In this study we investigate, via numer- 
ical simulations, the role of RT effects such as shadowing, self-shielding and filtering 
of the ionizing radiation, in giving raise to a fluctuating UVB. We focus on possible 
detectable signatures of these effects on quantities derived from Lya forest spectra, 
as photoionization rate fluctuations, r](= Ns e n/Nm) parameter distributions and the 
IGM temperature at z « 3. We find that RT induces fluctuations up to 60% in the 
UVB, which are tightly correlated to the density field. The UVB mean intensity is 
progressively suppressed toward higher densities and photon energies above 4 Ryd, 
due to the high He n opacity. Shielding of overdense regions (A ;> 5) from cosmic 
He ii ionizing radiation, produces a decreasing trend of rj with overdensity. Further- 
more we find that the mean rj value inferred from H i -He n Lya forest observations 
can be explained only by properly accounting for the actual IGM opacity. We outline 
and discuss several implications of our findings. 

Key words: cosmology: theory - radiative transfer - methods: numerical - cosmology: 
diffuse radiation - cosmology: large scale structure of Universe 



1 INTRODUCTION 

The lack of the Gunn-Peterson trough in the spectra of 
quasars at redshift z < 5 indicates that the intergalactic 
medium (IGM) after that epoch is photoionized by the ubiq- 
uitous presence of a metagalactic ionizing radiation, also 
known as the ultra-violet background (UVB). 

Usually, studies of the so-called Lya forest (the highly 
ionized absorbers that give rise to the absorption lines in 
quasar spectra) assume the UVB to be spatially uniform, 
with the same intensity and spectral shape everywhere in the 
universe at a given redshift. This is a good approximation if 
the sources of the UVB radiation are uniformly distribued in 
space, as long as the attenuation volume is large enough to 
contain a fair sample of both the sources and the attenuat- 
ing structures (Zuo 1992). This last condition would imply a 
low cosmic optical depth, r, to ionizing photons after reion- 
ization (z < Zi). The cosmic optical depth should be low 
enough infact, to insure the attenuation lengh to be larger 
than the typical separation between sources producing the 
UVB, which is \ a ~ 100/i -1 Mpc (comoving) for quasars 
and more than an order of magnitude less for galaxies. The 
relative contribution of these two populations to the UVB is 



unclear. Different independent studies found that the contri- 
bution from galaxies should be at least comparable to that 
of quasars (Giallongo etal. 1997; Shull etal. 1999; Steidel 
etal. 2001; Kim etal. 2001, Bianchi etal. 2001, Sokasian etal. 
2003) in the redshift range 2-3. As a caveat, we recall that 
Meiksin (2005) has pointed out that QSOs might dominate 
the UVB in the redshift range 2 < z < 3; according to his 
study QSO counts are highly uncertain near z = 3 and seem 
to be sufficient to account for the UVB at z = 2. The ex- 
pected signature of a significant contribution from galaxies 
is a considerable softening of the UVB spectrum compared 
to previous models assuming quasars as the sole sources of 
ionizing radiation (Haardt & Madau, 1996). 

Although the IGM at z < Zi is highly ionized, cosmic 
mildly overdense regions which give rise to the Lya forest 
contain enough neutral hydrogen and singly ionized helium 
to significantly attenuate the Lyman continuum flux from 
the sources of the UVB. As a result, the condition A p S> 
A s might not be fulfilled. For this reason, the metagalactic 
ionizing flux has been modeled starting from the intrinsic 
spectra of the ionizing sources, subsequently filtered through 
the IGM (Haardt & Madau, 1996; Fardal, Giroux & Shull, 
1998). In brief, these calculations consist in reprocessing the 
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UVB photons in a clumpy photoionized universe, by solving 
the radiative transfer (RT) equation in one dimension. The 
space density along the photon propagation direction, and 
physical properties of the sources (as well as of those of 
the absorbing clouds) are assumed as the average of the 
observed values, and the spectral redshift evolution of the 
UVB, assumed to be uniform in space, is derived. 

Clustering of the ionizing sources and scattering in their 
emission properties, together with inhomogeneities in the 
IGM density field can introduce significant fluctuations both 
in the intensity and spectral shape of the UVB. 

Indeed, increasingly strong observational evidence for 
a significantly variable metagalactic ionizing radiation 
has been accumulating through a number of studies 
of the H i and He n Lya forest. In a highly photoion- 
ized gas the ratio of He n to Hi column density, de- 
fined as the 77 parameter, is proportional to the ratio 
of H i to He n photoionization rates (Fardal etal. 1998): 
77 = NHeii/Nm oc r H i/r H eii- 

Kriss etal. (2001) and Shull etal. (2004) analyzed the 
FUSE observations of the fluctuating He 11 absorption to- 
ward the bright quasar HE 2347-4342 at z = 2.885. Us- 
ing H 1 Lya data from Keck and HIRES (High Resolution 
Echelle Spectrograph), they found an H 1 counterpart for 
more than 50% of the He 11 absorption lines. The measured 
77 values, cover a wide range from 1 to 1000 (with an aver- 
age (77) m 80). Analogous observations of the He 11 Lya for- 
est towards HS 1700+6416, by Remiers etal. (2004), reveal 
a similar 77 variation. The large scatter is a clear signature 
from an ionizing background which is significantly fluctuat- 
ing throughout the IGM. Shull etal. (2004) pointed out two 
important points: (i) a small scale 77 variability (on typical 
scales of Az ~ 10~ 3 , corresponding to about lh" 1 Mpc co- 
moving at z ~ 3), and (ii) an observed correlation of high-77 
(i.e. soft ionizing spectra) absorbers with low density regions 
(voids in the H 1 Lya distribution) . The authors further sug- 
gest that these effects, confirmed by Reimers etal. (2004) ob- 
servations, might be caused either by spatial/spectral fluc- 
tuations of the ionizing sources on scales of 1 Mpc, or by RT 
effects trough a filamentary IGM whose opacity variations 
control the penetration of 1-5 Ryd radiation over 30-40 Mpc 
distances (a combination of the two is also possible). 

Although recent HST (Hubble Space Telescope) obser- 
vations (Telfer etal. 2002) show a broad distribution of QSO 
spectral indices from a ~ to 3 (J v oc v~ a representing the 
quasar spectrum) this spread could hardly be advocated to 
explain the small scale 77 fluctuations. The reason is that 
the QSO mean separation (« 100 Mpc comoving at 2 ~ 3) 
exceeds by far the characteristic scale of about l/z _1 Mpc 
on which the rj variations are observed. Radiative transfer 
through an inhomogeneous medium appears to be the most 
natural origin for such fluctuations. In turn, these might 
result from a number of effects, such as shadowing, self- 
shielding and filtering of the radiation. All these effects act 
on the same scales of the inhomogeneities in the IGM, thus 
on scales smaller than 1 Mpc. It is then important to study 
the correlation function of both J v and 77 with the density 
fluctuations. 

The first attempt to account for the RT effects of a 
background ionizing radiation was by Nakamoto, Umemura 
& Susa (2001). They focused on the reionization of an in- 
homogeneous universe looking at the evolving configuration 



of the neutral hydrogen distribution, but not directly at the 
fluctuations in the radiation field. Nonetheless, emphasizing 
the deviations of the neutral hydrogen distribution obtained 
in their simulation from the case of a perfectly transparent 
medium embedded in a uniform ionizing background, they 
indirectly stressed the importance of RT effects such as self- 
shielding and shadowing by translucent regions, which in- 
deed have been well reproduced by their simulations. 

Radiative transfer is furthermore expected to have some 
effect on the IGM temperature. Previous studies (Abel & 
Haehnelt 1999; Bolton etal. 2004) focused on the effects of 
RT on the temperature evolution during reionization, i.e. 
behind an ionization-front (I-front) expanding in a neutral 
or mildly ionized IGM. These analysis have shown that ra- 
diative transfer can boost the mean gas temperature with 
respect to the case in which an instantaneous, uniform reion- 
ization occurs (i.e. its OT counterpart), as a consequence of 
filtering effects. In fact, the ionizing spectrum across the 
I-front becomes harder as it expands further from the ioniz- 
ing source, because most of the low-energy photons (having 
higher cross-sections) are absorbed by recombining atoms 
inside the H 11 region, and its photo-heating power increases 
accordingly. The overall predicted RT effects are: (i) an en- 
hancement of the IGM temperature by a factor ~ 2 after 
reionization; (ii) increased scattering in temperature-density 
relation; (iii) a strong hardening of the metagalactic ioniz- 
ing flux after the He 11 reionization. Here we consider the 
more stationary situation in which the mean ionization level 
in the simulation volume remains roughly constant; hence 
strong RT signatures on the IGM temperature are not ex- 
pected. Nevertheless we look for a systematic dependence of 
the photoheating rates on the density. 

So far, much of the work in studing UVB fluctuatios 
have be done to asses the effects produced by inhomogeneous 
distribution of ionizing sources, almost neglecting or just 
approximately including the effects of RT through the fila- 
mentary IGM. Zuo (1992a,1992b) first developed a theory to 
deal with the fluctuations produced by randomly distributed 
ionizing sources, in terms of the probability distribution of 
the ionizing UV radiation field (Zuo, 1992a) and of its two- 
point correlation function (Zuo, 1992b). Meiksin & White 
(2003) made a step forward by appying Zuo's formalism to 
ionizing sources distribution derived via numerical simula- 
tions of structures formation. A few groups have studied, 
trough numerical simulations, the effects on the statostical 
properties of the Lya forest of UVB fluctuations induced 
by source clustering or non-uniform emissivity properties as 
light-cone effects, QSO time variability, lifetime and spec- 
tral index scattering. Gnedin & Hamilton (2002) first inves- 
tigated via RT simulations the effects of UVB fluctuations 
induced by quasars inhomogeneous distribution, on the mat- 
ter power-spectrum inferred from the observed/simulated 
trasmissivity of the Lya forest (Croft etal. 1999). By com- 
paring the mock spectra from RT simulations and from sim- 
ulations assuming photoionization equilibrium with a uni- 
form UVB, they found the effects of RT to be negligible for 
the power spectrum; nevertheless it is interesting to notice 
that, in order to match the observed transparency of the 
IGM, the mean photoionization rate required when includ- 
ing properly RT, has to exceed by 20% the one required in 
the OT approximation. Meiksin & White (2004) attempted 
to quantify the UVB fluctuations and examined the possible 
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effects on the power spectrum and auto-correlation function 
of the Lya forest. They calculate the UVB intensity smear- 
ing out the radiation emitted from a discrete distribution 
of QSOs in the simulated volume, within an attenuation 
length whose evolution with redshift is evaluated modeling 
the gas density distribution with a PM simulations. Their 
calculation separately adds the contribution of Lyman Limit 
System (LLS), which are not resolved in numerical simula- 
tions, on the basis of their observed statistical properties. 
Croft (2004) uses a time-independent ray tracing algorithm 
to model the space density of ionizing radiation produced 
by QSOs at z » 3, neglecting its feedback on the ioniza- 
tion state of the IGM. The study is based on N-body sim- 
ulations and prescriptions must be used to relate the dark 
matter density and the H i optical depth. McDonald etal. 
(2004) use a rough self-shielding approximation, in which 
the background radiation in each cell is attenuated by the 
column density of that cell; shadowing and filtering effects 
are neglected altogether. The most relevant outcome of these 
studies for the present one is that UVB fluctuations induced 
by the discreteness of the UVB sources could be only rele- 
vant at high redshifts and on large scales. Below z = 4 the 
importance of such fluctuations rapidly decreases, as a con- 
sequence of the increase of A p which smooths the flux at a 
given position over a larger ensemble of sources. 

For this reasons we take here a complementary, and still 
unexplored, approach which consists in the investigation of 
UVB fluctuations introduced by the RT effects mentioned 
above. We focus on possible detectable signatures of these 
effects on quantities derived from Lya forest spectra, as pho- 
toionization rate fluctuations, the temperature-density rela- 
tion and f] parameter distributions. To this aim we have 
performed RT simulations, applying CRASH (Maselli, Fer- 
rara & Ciardi, 2003, hereafter MFC03) to a cosmological 
density field exposed to an ionizing background radiation 
characterized by a spatially constant emissivity. Such radi- 
ation field, by interacting with the filaments of the cosmic 
web, is rapidly turned into a fluctuating one. The rest of the 
paper is devoted to study and quantify this effect. 



2 NUMERICAL SIMULATIONS 

The aim of the present study is to estimate the effects of 
RT on the metagalactic ionizing radiation and the resulting 
effects induced in the physical properties of the IGM. The 
present study is focused at z w 3. We have exposed a snap- 
shot at a fixed redshift of a cosmological density field to 
an (initially) uniform ionizing radiation field. The density 
field has been derived with the numerical code described 
in Marri & White (2002), a modified, multiphase version of 
the entropy-conserving SPH code GADGET2 (Springel etal. 
2001). Several algorithms have been designed to reduce ar- 
tifacts which occur in cold, dense gas clouds embedded in 
a hot diffuse halo; the code includes a new implementation 
of feedback which allows supernova energy to be channeled 
effectively into the heating of diffuse gas and the evapora- 
tion of cold clouds. A multiphase description of the inter- 
stellar medium is adopted, based on an explicit separation 
of protogalactic gas into diffuse and dense (star forming) 
components which is able to suppress star formation, reheat 
cold gas and drive outflows from galactic disks. The pres- 



ence of a UV background produced by QSOs is included in 
the OT approximation following Haardt & Madau (1996) 
from which the IGM temperature and ionization state are 
calculated assuming photoionization equilibrium. 

We use here a snapshot at z = 3.27 from a simulation 
which has been used and analyzed in previous works of our 
group (Bruscoli etal. 2003; Maselli etal. 2003; Fangano, Fer- 
rara & Richter 2005, in preparation). These studies have 
shown that the statistical properties of the IGM, e.g. the 
PDF of the transmitted Lya flux, are correctly reproduced. 
The simulation uses 128 3 particles in a 10.5ft _1 comoving 
Mpc box (corresponding to a mass resolution of 8.5 x 1O 6 M 
in gas) and assumes a ACDM cosmological model with 
Q = 0.3, Q A = 0.7, Q b = 0.04, h = 0.7 and cr 8 = 0.9; 
periodic boundary conditions are assumed. 

2.1 Radiative Transfer 

We have performed the RT calculation using the code CRASH 
. A detailed description of the code, its performances and 
several test problems can be found in MFC03, to which we 
defer the interested reader. In the present study we neglect 
the effects of cosmological expantion; this choise is moti- 
vated by the fact that the light crossing time of an attenu- 
ation leght (~ 100 Mpc) is roughly one tenth of the Hubble 
time at the redshift of interest. We have exposed the density 
field obtained from the SPH simulation to a UVB radiation 
with spatially constant emissivity. The UVB is subsequently 
modified as it propagates in the filamentary density field by 
self-shielding, shadowing and spectral filtering effects, finally 
resulting in a fluctuating field. The RT simulation has been 
evolved for a physical time, t s = 0.7 x 10 7 yr. This choice 
has been made, after some experiments, to guarantee that 
numerical convergence of the solution is achieved. It has to 
be noticed that here the simulation time does not corre- 
spond to a physical evolution of the background, but it is 
the time required for the gas ionization and temperature 
to reach equilibrium in the self-consistently radiation field 
shaped by radiative transfer. In the following we discuss two 
crucial, albeit somewhat technical, points concerning a new 
implementation of a uniform radiation background in CRASH 
and its flux calibration within the simulation. 

2.1.1 Simulating a uniform UVB 

The algorithm described in MFC03 for the implementation 
of the presence of a uniform background radiation in the 
simulation box has been improved in several ways. In the 
original implementation, background radiation was modeled 
by allowing box boundary cells to emit photons packets with 
the prescribed UVB spectral energy distribution. Although 
this choice is suitable for some type of problems it cannot 
match the high accuracy required by the present study. To 
better fulfill such requirement we have devised a new im- 
plementation in which all cells in the box, but the ones 
with overdensity A — p/{p) ^ A e , behave as UVB radia- 
tion sources, in the sense that they emit photon packets 1 . 
The threshold on the overdensity is meant to preserve the 

1 It is important to notice that here the cells emitting photons 
are not meant to reproduce physical sources of UV radiation; 
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effect of self-shielding in high density regions and galactic 
halos. For definiteness, we take the value corresponding to 
the typical overdensity of collapsed halos at the virial ra- 
dius. The exact value of A e depends on the density profile 
of the halo; it can be shown that it is w 59.2 for a isother- 
mal profile and it is » 63.7 for the NFW profile (Navarro, 
Frenk & White 1996). In this paper, A e = 60 has been 
assumed. We adopt periodic boundary conditions, so that 
packets can travel trough the simulation box for a number of 
cycles, N cyc , given by the ratio between the mean free path 
of ionizing photons, A p , and the physical size, L, of the sim- 
ulated volume. We have derived an estimate for X p , running 
a benchmark simulation aimed at evaluating the statistical 
distribution of the mean free paths of photons -whose fre- 
quencies are distributed according to the HM96 spectrum- 
in the assigned density field. For this purpose, in the bench- 
mark simulation we just follow the propagation of photons 
packets, i.e. not cosidering their effect on the gas properties; 
due to the statistical treatment of the packet energy depo- 
sition inherent to the Monte Carlo scheme (see MFC03), 
we define as "absorbed" a packet whose initial energy has 
been depleted by > 95%. In Fig. 0we show the mean free 
path distribution obtained from such numerical experiment. 
The solid curve corresponds to the mean free path obtained 
by including the opacity of both hydrogen and helium ab- 
sorbers, whereas the dotted curve corresponds to H i opacity 
only. The mean (median) values of the two distributions are 
300 and 495 (199 and 306) h' 1 Mpc comoving. This shows 
that the He n contribution to the diffuse IGM opacity is 
significant. The above value of X p , for H i opacity, is consis- 
tent with that inferred from Fan etal. (2002), X p « lOO/i" 1 
Mpc at z ~ 3, by comparing observations of high resolution 
quasar spectra with semianalytical models of the Lya forest 
with a simple schematic model of the IGM opacity based on 
a peak statistics approach. 

Guided by this finding, and given our box size, we fix 



N c 



10. The agreement of the simulated and observed 



mean values of X p further supports the N cyc choice above. 
Periodic boundary conditions allow us to simulate volumes 
smaller than the mean free path of neutral hydrogen ioniz- 
ing photons, having thus a higher spatial resolution, with- 
out much information loss. However, we make clear that the 
main aim of the present study are RT effects on scales of 
~ 1 Mpc, the characteristic scale of IGM filaments; no at- 
tempt is made to describe large scale fluctuations induced 
by non-uniform distribution of the UVB sources. 

According to the CRASH Monte Carlo algorithm track- 
ing photon packet propagation, the mean number of cross- 
ings events per cell is given by 



N cyc N p 
Ni 



(1) 



where N p is the total number of photon packets in a sin- 
gle run, and iV 3 the number of grid cells in the computa- 
tional volume 2 . The new background implementation has 
been tested by running a simulation in which we force the 



the adopted method aims at reproducing a uniform UVB in the 
simulated volume, within a Monte Carlo ray tracing scheme. 
2 In the OT case, a single packet crosses typically N cyc N c cells, 
thus the total number of crossings in a simulation is N p N cyc N c 
for JVf cells (for more details see MFC03) 
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Figure 1. Distribution of the comoving mean free path of ionizing 
photons at z = 3.27 as derived from the SPH simulation. The solid 
curve includes opacity of both H I and He II ; dotted curve is for 
H I only. 



condition r = 0. Although, when not inlcuding cosmological 
expansion, the r = assumption would imply an infinite 
radiation filed (i.e. the Olber's paradox), we get around this 
problem by fixing N cyc to a finite value. All the parameters 
then, are setted in order to reproduce the required UVB in- 
tensity, as described in Sec. 2.1.1. Ideally, in this case the 
intensity of the UVB should be perfectly uniform, as fluctua- 
tions caused by photon-matter interactions (absorption and 
scattering), redistributing the ionizing radiation in space 
and frequency, are suppressed. In practice, uniformity can 
be achieved only up to a certain accuracy level, due to finite 
number packet statistics, which produces noise fluctuations 
in the emission pattern. The amplitude of such fluctuations 
decreases with increasing N p , which is however constrained 
by computational expense. We have performed different test 
runs increasing N p to accurately quantify the amplitude of 
such fluctuations. The computed r.m.s. of the fluctuations, 
<Tjv, around the mean value of the UVB intensity as a func- 
tion of the numerical parameters of the run are shown in 
Table 1. The value for N c and N cyc are the same adopted 
for the actual cosmological simulations 3 . These fluctuations 
are due to the photon shot noise intrinsic to the adopted 
numerical algorithm. 

From Tab. 1, we see that runs at higher resolution (i.e. 
higher values of N p ) are characterized by a small intrinsic 
scatter of the radiation intensity: ajv, decreases from 20% 
for N p = 3 x 10 5 , to less than 1% for N p = 3 x 10 8 . At 
the same time, we have also checked that during these test 
runs (Afcross) closely matches the theoretical value (as ex- 
pressed in eq. . As a good compromise between accuracy 
and computational time we will adopt N p = 3 x 10 7 for the 
cosmological RT simulations. This value, according to the 
above tests, guarantees that the ionizing radiation field fluc- 
tuations are uniformly distributed within a precision level of 

3 The intrinsic scatter crjv , associated to the numerical resolution 
of a given run, does not depend on the physical time of the sim- 
ulation, t s , and on the intensity of the UVB which, as discussed 
later, depends on the normalization factor A. 
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Table 1. Parameters of the simulation test runs used to check uniformity and calibration of the radiation field. 
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Figure 2. Spatial fluctuations of H I (left panel) and He II (right) photoionization rates in a slice across the RT simulation of (physical) 
depth Ax fa 27 kpc. The curves denote the mean values. 



4%. Hence, larger fluctuations seen in the simulation outputs 
can be safely attributed to RT effects. 



2.1.2 UVB intensity calibration 

Particular attention must be devoted to the ionizing flux 
calibration. The problem arises when the emissivity has 
to be related to the photon content of packets originat- 
ing from a cell, as it is not obvious how convert the latter 
into the specific flux units in which J„ is usually assigned 
(erg s _1 cm" 2 Hz" 1 sr" 1 ). To overcome the problem we 
have devised a calibration technique which turned out to be 
highly precise and straightforward. 

Given a set of simulation parameters, N p , N cyc , N c and 
t a , the intensity of the radiation field can be tuned by vary- 
ing the number of (monochromatic) photons contained in a 
packet, jV 7 , at emission. We express this quantity as follows 



iV 7 = A At 



Jiy) 
J912 



= A 



N*ts J{v) 

Np N C yc J912 



(2) 



where A is the calibration factor 4 to be determined and At 
is the mean time interval between two subsequent cross- 
ings of a cell by photon packets; for the second equality see 
MFC03. A nice feature of the previous expression is its scal- 
ing with the simulation parameters. Once the intensity and 
the UVB spectral shape are fixed, the factor A is the same 
for any combination of the numerical parameters. An order 
of magnitude for A can be derived by equating the theoreti- 
cal and numerical estimates of the ionization fraction incre- 
ment, Axmi, for a cell during the time interval At. We will 
consider here, for the sake of clarity, the case of a pure hy- 



4 Note that the calibration factor A has physical dimensions s 1 , 
i.e. a photoionization rate. 
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drogen gas and a monochromatic spectrum (J(v)/ J912 = 1). 
The theoretical value is given by 

Axhii = TrnxmAt; (3) 

similarly, the numerical one can be expressed as 

Aa&i = JV 7 (1 - e~ T ")/N m = AAtr c /N m (4) 

where r c is the cell optical depth, Nm = iihiAi 3 , with Air 3 
being the cell volume, is the total number of H-atoms in the 
cell; the last equality holds for the OT case (r c <C 1). By 
equating the two expressions, one can write A as 

A = ^Ax 2 , (5) 

where we have used r c « nmooAir, with ao being the 
H 1 photo-ionization cross-section at the Lyman limit. For 
the RT simulations used in this study, Air ~ 27 kpc; hence 
for a typical r m = KT 12 s _1 , A » 10 51 s" 1 . 

The determination of the exact value of A for the gen- 
eral case of a H/He gas and of a non- monochromatic spec- 
trum, is not as straightforward. The adopted strategy has 
been to select a cell at the center of the simulation box, to 
assign it physical properties characteristic of the hydrody- 
namic simulation (we choose n = 10 -5 cm -3 , T = 10 4 K, xhi 
the value at the photo-ionization equilibrium with the as- 
signed Thi value). We assume r = for all remaining cells 
in the box; this guarantees that the ionizing radiation field 
impinging on the target cell has not been modified either 
in intensity or spectral shape by RT effects. In the above 
configuration, we run CRASH to compute the H 1 fraction in 
the central cell and the average photoionization rate, Thi, 
by counting the cumulative number of ionizing photons ab- 
sorbed in a given simulation time. We then tune the value 
of A by requiring that the Xbi and Fhi estimated values 
match the assigned values within an accuracy of 1%. We 
find that r H i = 1.2 x KT 12 s" 1 (HM96) corresponds to 
A = 1.55 x 10 52 s _1 . 

As discussed above, the same value of A reproduces cor- 
rectly the assigned ionization rate (for a given UVB spectral 
shape) independently of the values of N p , N cyc , N c and t 3 , 
provided that 7V 7 is evaluated according to eq. [5] The de- 
gree of accuracy is shown in the bottom part of Tab 1 where 
the r.m.s. deviation, crjv, of the distribution of Aa;g n with 
respect to the theoretical value Aa;jjn i s shown for different 
combinations of the numerical parameters. The scatter is al- 
ways less than 10%, depending on the resolution of the run. 
The calibration method adopted results thus to be accurate 
and stable. 

In actual RT simulations, we should further account for the 
opacity of the IGM in order to reproduce at each cell the 
correct ionizing flux 5 ; the photon content of packets is infact 
depleted while it propagates trough the box, due to non-zero 
opacity of the IGM distribution. To counteract this effect, 
we boost the normalization parameter A by 2.3; the exact 
value of this factor has been determined empirically, requir- 
ing that the mean opacity of the simulated volume is the 
same for the OT (zero-opacity) and RT cases. If the calibra- 
tion factors for the RT and OT cases are related according 



5 With "correct" we mean the ionizing flux which reproduces the 
assigned mean opacity. 



to the equation A = A e~ T ' f f , it corresponds to an 
effective optical depth r e a ~ 0.8. 

Extracting the correct physical values of the photoion- 
ization and photoheating rates (at each cell) from the RT 
simulation is not straightforward as well. Again, we need a 
reference OT simulation (with the same numerical param- 
eters) in order to calibrate correctly these values. The OT 
and RT simulations differ (i) in the normalization factor, 
A OT and A RT related as above, and (ii) in the fact that for 
the OT we do not allow for absorption to occur. In both RT 
and OT cases, summing up the contribution given at each 
crossing event, the following relations hold: 

^ = a N -r(n)°i(n)> ( 6 ) 

j Across 

Hi = C'i Nj(u-i)ai(uj)(uj - Vth.i), 

j Across 

where j is the index counting the crossing events in a given 
cell and Uj is the frequency of the packet at the j'-th crossing; 
d and C[ are a scaling factor which depend only by the 
numerical parameters and are the same for both the OT 
and RT cases. The other symbols follow the conventional 
notation. Having determined A OT as seen above, we can 
safely impose that (if T ) = Yf {e.g. rg T = 1.2 x lO'^s -1 ) 
and accordingly derive C, which in turns is used to derive 
the correct Tf T physical values at each cell in the box. 

Finally, we mention that in the RT simulation the emis- 
sion of diffuse radiation from recombining gas (estimated to 
account for about 30% of the ionizing background at z = 3) 
has been turned off inside CRASH. This is done to prevent 
double-counting of such radiation, which is already included 
in the Haardt & Madau (1996) computation of the adopted 
UVB. 



3 RESULTS 

In this section we present the results of the RT simulations. 
We first analyze the amplitude of RT-induced fluctuations as 
traced by UVB photo-ionization and photoheating rates of 
the three absorbers H 1 , He 1 and He 11 . These are defined 

as 

5ri =-(f^-' 5Hi =— wr' (7) 

the index i £ {H 1 , He 1 , He 11 } denotes the absorbing 
species and the averages are done over all the cells at the 
end of the RT simulation We then focus on the effects of 
such fluctuations on the IGM r\ = Nnell/Nm. distribution, 
the ^-density relation and on the IGM temperature-density 
relation, with the aim of identifying possible detectable sig- 
natures. 

3.1 UVB fluctuations 

Spatial fluctuations of the H 1 and He n photoionization 
rates are shown in Fig. [5] on a slice across the RT simu- 
lation box of (physical) depth Air ~ 27 kpc, corresponding 
to one grid cell (<5Fhi and 5Fh c ii in the left and right pan- 
els respectively). RT effects result in significantly fluctuat- 
ing photo-ionization rates on scales of few comoving Mpc, 
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Figure 3. Spatial fluctuations of H I (left panel) and He II (right) photoionization rates, as a function of the local overdensity, A. Each 
point in the figures corresponds to a cell of the simulated grid; the solid lines show the mean SV value as a function of A. 
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Figure 4. Probability distribution function of the fluctuations for three selected values of the overdensity, A = 0.3 (solid line), A = 1 
(dashed) and A = 3 (dotted-dashed). The left and right panels refer to <5rjji(A) and <5rjj c ii(A), respectively. The dotted curves are the 
analogous distributions at A = 0.3 for the OT case. 



whose structure is closely linked to the underlying density 
fluctuations. As shown in the following, typically regions of 
low overdensity (say A ^ 1, often referred to as voids) have 
5T > 0, whereas overdense regions tend to have negative 
fluctuations. Comparing the two maps one can also notice 
that ThcH fluctuations are much stronger than those of Thi- 
In general, though, RT induces fluctuations of several tens of 
percent of the mean ionizing flux in both species. Such devi- 
ations from a homogeneous radiation field are produced by 
the combined action of shielding and shadowing effects. Al- 
though these features can be qualitatively understood from 
simple physical arguments, their quantitative evaluation can 
only be done through full 3D RT simulations. 
In order to make the above statements more quantita- 
tive, we show in Fig. [3] the <5Fhi and <5rHeii distributions 
as a function of the local overdensity, A. Each point in 
the plots corresponds to a cell of the simulated grid. The 
plot in left panel shows clearly how the mean amplitude 
of H i photoionization rate fluctuations varies with overden- 



sity, passing from positive values (of the order of a few per- 
cent) in the voids to negative ones, larger than 10%, towards 
higher densities. <5rHeii, shown in the right panel, exhibits 
an analogous trend, although with larger excursion, with 
values ~ +25% (~ —50%) at the lowest (highest) density 
boundaries of the density range. Deviations from the mean 
value can be as high as 20% and 60%, for Thi and ThoII re- 
spectively, in structures corresponding to overdensities with 
A J> 50, where therefore the UVB intensity appears to be 
significantly suppressed. The steeper dependence and larger 
scatter of ^ThoII-A relation reflect the fact that the uni- 
verse is more opaque in He n than in H i . Note also that in 
underdense regions FhcII can exceed by more than 40% the 
volume averaged mean value, (THen). 

We focus now on mildly overdense regions, such as filaments, 
which give raise to absorption lines in the Lya forest. These 
regions are exposed to roughly the mean photoionization 
rates for H i , but the photoionization of He n is here sup- 
pressed by about 20% with respect to the mean value. This 
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fact suggests that the He in fraction, which both in numeri- 
cal simulations and data interpretation of Lya forest is usu- 
ally derived assuming photoionization equilibrium with a 
uniform UVB , is probably seriously overestimated. The ul- 
timate physical nature of the decreasing trend of 8Ti, with 
i £ {H i ,He n }, with density is a mixture of self-shielding 
and shadowing effects. 

It is very instructive to inspect the probability distribu- 
tion function of the fluctuations at fixed overdensities. Fig. 
2] shows such distribution at A 6 {0.3,1,3} (solid, dashed 
and dotted-dashed curves respectively), for 5Fhi(A) (left 
panel) and <5FHeii(A) (right panel). The distributions re- 
semble Gaussian curves whose mean value decreases with 
A. The shift of the mean is caused by the increasing impor- 
tance of the shielding and it is about five times larger for 
He ii (= —0.25) than for H i (= —0.05): as explained above, 
the attenuation of the flux due to the structures shielding 
is more prominent above 4 Ryd due to the larger opacity at 
those energies. In the case of H i , a significant increase of 
the distribution width is observed towards larger overdensi- 
tites; contrarily, the three distributions obtained for (SFhcII 
show a similar dispersion. This results from the fact that, 
while the H i optical depth of the IGM is dominated by rare, 
overdense regions, the He n optical depth is non-negligible 
even in voids. This fact is confirmed if one look at the proba- 
bility distribution for the mean free path of photons at 1 Ryd 
and 4 Ryd, shown in Fig. |S] (solid and dotted-dashed curves 
respectively). For 1 Ryd photons the mean free path distri- 
bution is highly peaked at roug hly 60 h' 1 Mpc 6 , showing 
that these photons are preferentially absorbed in high den- 
sity peaks. The 4 Ryd photons mean free path distribution 
is wider and has a mean value of about 7 ft -1 Mpc, which 
is a clear indication for a significant contribution of voids to 
the He n opacity. 

We have checked that the dispersion of the flux 
distribution shown in Fig. 2] do not result from spurious 
effects introduced by the emission algorithm, comparing the 
distributions above with the analogous ones derived in the 
OT simulation 7 . The latter are shown in the panels of Fig. 
|1] by the dotted curves for A = 0.3 (as the OT distribution 
does not depend on overdensity, we only plot this case). For 
visualization purposes, these curves have been normalized 
so that their peaks match the peaks of the corresponding 
RT distributions. The dispersion of these distributions is 
obviously unphysical as it represents the intrinsic error of 
the numerical algorithm, which is of the the order of 4% 
for H i , as already pointed out in Table 1, and more than 
twice for He n (~ 10%); the higher numerical dispersion for 
He ii is due to the worse spectral sampling of energies above 
4 Ryd 8 . Both for H i and He n , the dispersion induced by 
RT in the photoionization rates values is significantly larger 
than the one originating purely from numerical noise and, 

6 Note that these values are lower than the ones in Fig.0 because 
here we consider only photons at the H I and He II ionization 
threshold. 

7 We use the OT simulation we have run to test the accuracy 
of the uniformity of the emission algorithm, with N p = 3 X 10 7 , 
having thus the same numerical resolution of the RT simulation. 

8 The number of photon packets emitted with energy above the 
He II threshold is roughly one order of magnitude smaller less 
than N p . 



in the case of H i , it is increasingly so at larger values of 
A. Hence, this confirm that the effects we are finding are 
of physical origin and the corresponding fluctuations are a 
truly genuine by-product of UVB filtering through cosmic 
density structures. 

We have found so far that RT effects induce significant 
fluctuations in the metagalactic ionizing radiation; these 
could be effective in deviating the actual IGM ionization 
state from the one inferred assuming a uniform ubiquitous 
UVB and hence should be accounted for in high-precision 
studies of the Lya forest. However, performing full RT simu- 
lations is extremely expensive and prohibitive in maximum- 
likelihood analysis. The tight correlation between fluctua- 
tions in the photo-ionization rate and overdensity can be 
used to correct the commonly adopted mean values derived 
in calculations (e.g. Haardt & Madau 1996; Fardal etal. 
1998) carried out with ID RT calculations, under the as- 
sumption of photo-ionization equilibrium with a uniform 
UVB. 

In order to cast our results in a form which could be useful 
for future studies, we have computed the photo-ionization 
rate correction factors from our RT simulations. These are 
defined as 

rf T (A) =/r 4 (A) T? T (A); (8) 

their mean values are plotted in Fig. HJ for H i , He i and 
He ii (upper, central and bottom panel respectively) along 
with the associated la errors. These results are valid for 
z w 3. In a companion paper we will study the evolution 
with redshift of the UVB fluctuations. RT transfer effects 
are generally more important in overdense regions, although 
the correction factors derived in the present study can be 
reliably used for A < 100, a limit dictated both by mass 
resolution of the SPH simulation and shot noise associated 
to the limited number of photon packets. Differently from 
the case of H i , for which (/r H i) ~ 0.98 is close to unity, the 
mean value of the correction factors for He i and He n are 
(/r HeI ) « 1.56 and (/r HeII ) ~ 0.2 (averages are evaluated 
over the overdensity range 0.03 < A < 100) . This reflects 
mean opacities differences among the three absorbers. In the 
OT approximation, the IGM opacity is higher in H i than in 
He i and, as a consequence, the UVB increase necessary to 
match the IGM H i Lya opacity in RT simulations results 
in a higher mean ThcI- For the same reason, as the IGM 
OT opacity is by far higher in He n than in H i , (FhoIi) 
is highly suppressed with respect to the OT value when 
RT effects are properly treated. Again, our results demon- 
strate that the contribution to the He n opacity of the dif- 
fuse IGM (as opposed to discrete absorbing systems) is rel- 
evant. The same arguments explain the steeper dependence 
of the He n correction factor from density and its larger dis- 
persion. 

A few remarks are worth concerning the correction fac- 
tors. First, we remind that the adopted flux calibration (see 
Sec. 2) guarantees that the mean H i opacity in the RT sim- 
ulation box corresponds to that at the photoionization equi- 
librium with the uniform UVB in the OT case. We found 
that the mean value of Thi in the two cases (RT and OT) 
are the same within few percent. Gnedin & Hamilton (2002) 
found a different result: in order to reproduce the same 
H i Lya IGM trasmissivity in their RT calculation, they 
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Figure 5. Probability distribution of the mean free path of pho- 
tons at 1 Ryd (solid curve) and 4 Ryd (dotted-dashed curve). 



require a 20% higher volume-averaged H i photoionization 
rate with respect to the OT case. The reasons for the dis- 
crepancy are not clear, but it is probably due to the fact that 
while Gnedin & Hamilton (2002) account for RT during the 
reionization process (i.e. ionizing a neutral gas), here we re- 
distribute the ionizing radiation accounting for RT effects in 
a pre-ionized gas. Meiksin & White (2003) also found that 
fluctuations required a boost in the radiation field, com- 
parable to Gnedin & Hamilton; however, in a more recent 
investigation, they found that the boost is reduced to few 
percents at z ~ 4 when correlations in the radiation field due 
to RT were included (Meiksin & White 2004). This seems 
to suggest that the boost depends on the magnitude of the 
correlations in the radiation field. 

It is furthermore interesting to note that Thi can be fit- 
ted with a power-law of the form: F H i(A) = r HI A~ ai 
with oi — 0.022. In the calculation of optical depth to the 
Lya scattering given by (e.g. Weinberg etal. 1999; McDon- 
ald & Miralda-Escude 2001) 



T 



(i + zf (n b h 2 ) 2 

T^H(z)Y- 1 2,m(^Y 



0) 



it is thus possible to account for UVB spatial fluctuations 
adopting a corrected index — /3+ai, i.e. r fx A^ /r^ 12 hi> 
i.e. in a variation of the apparent IGM equation of state. 
Even though at z fx 3 we found a negligible variation of j3 
(only few percents), we expect the above effect to be impor- 
tant at higher redshifts. 

Fig-Qshows the PDF of the hydrogen photo-ionization 
rates resulting from our simulations (normalized to the vol- 
ume averaged mean value (Fhi}): the solid curve refers to 
the full RT simulation, whereas the dashed one has been ob- 
tained from the OT simulation. The RT distribution is much 
wider than the OT one, showing an extended tail towards 
lower values of Fhi /(Fhi) which is caused by shadowing and 
shielding due to gaseous structures. There is also a slight 
increase of cells seeing an ionizing flux above the mean, oc- 
curring as a result of the flux amplitude redistribution. Fur- 
ther interesting aspects can be inferred comparing the UVB 
fluctuations produced by RT effects and those induced by 
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the inhomogeneous distribution of rare quasars, assumed to 
be the dominant UVB sources, derived by Croft (2004) and 
shown in Fig. 5 of that paper. Even in that case, the PDF 
of the ionizing flux appears to be asymmetric: however, the 
asymmetry is in the opposite direction. A similar result is 
obtained by Meiksin & White (2003) as shown in Fig. 4 of 
that paper, where the effects of the sources stochastic dis- 
tribution in space are studied. In fact, fluctuations origi- 
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nating from sources inhomogeneous distribution are charac- 
terized by considerable power in the high flux tail (j> 1.5 
times the mean value), corresponding to the enhanced ion- 
izing flux in the vicinity of the luminous ionizing sources. 
Hence, both RT and inhomogeneous distribution of ionizing 
sources, contribute to broaden the radiation field distribu- 
tion, basically increasing the power on opposite sides of the 
flux mean value. It is important to realize though, that the 
spatial character of these fluctuations is very different. The 
signature of source inhomogeneities is seen on scales compa- 
rable to the typical separation of the sources, i.e. > 100 Mpc 
for quasars. RT-induced fluctuations are instead found on a 
few (comoving) Mpc scales. To clarify this point we have 
studied the spatial correlation of RT-induced UVB fluctua- 
tions, computing the SFi correlation function by using the 
following formula: 



= 1 k^j 



SFjjj) SFjjk) 



(10) 



Here the index j and k span all the cells in the simulation 
grid, A is the comoving separation between the j-th and k- 
th cells and A^a the number of pair cells whose distance is 
the range [A, A + AA] . We have similarly calculated a second 
correlation function, C[(A), obtained from a random realiza- 
tion of SFi, in the same volume and with the same mean as 
Ci(A). The difference, C<(A) = Ci(A)-Cf(A) is plotted in Fig. 
[S]for H i and He n . The zero-point of these functions gives 
the characteristic correlation comoving length of the RT in- 
duced fluctuations, which is found to be roughly 5.5 and 
3.5 Mpc comoving for H i and He n respectively. Although 
this results should be checked for convergency, by perform- 
ing analogous simulations in larger volumes with the same 
resolution (Meiksin & White 2004), they seem to be plausi- 
ble. The correlation scales we have found, could be associ- 
ated with the typical size/separation of filaments (i.e. mildly 
overdense regions) of the IGM which are clearly identified in 
the maps presented in Fig. These structures have optical 
depths > 0.1 both in hydrogen and helium, and therefore 
can leave their imprint in the fluctuation field. 



Meiksin & White (2003,2004) have shown that inhomoge- 
neous sources (assumed to be QSOs) distribution can pro- 
duce correlations which are significant even on small scales. 
Both the RT induced fluctuations and the ones originat- 
ing from disomogeneous sources distribution, should then 
be treated at once, in order to derive a correlation function 
which is reliable for the study of the pixel fluxes statistic in 
high-z QSOs spectra. 

3.2 He ii /H i ratio fluctuations 

The rj parameter, defined as the ratio of He n to H i column 
density of an absorber in the Lya forest, gives informations 
on the intensity and on the spectral properties of the meta- 
galactic ionizing radiation. Assuming photoionization equi- 
librium in highly ionized H/He gas, the following expression 
can be derived (Fardal etal. 1998): 
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rate coefficients, appropriate in the low density regime. We 
have derived the rj parameter, for each cell of the numeri- 
cal grid, from eg. 1111 with WHeHi/rami and rHi/FHeii taken 
from the RT simulation. The results are plotted in Fig. |5] 
the solid (dotted) line represents the mean rj value derived 
in the RT (OT) simulation and the dashed curves are the 
lcr r.m.s. deviations of the distribution. 

In the OT case, as expected , rj shows little dependence 
on A as a result of the fact that the three ratios in eq. 1111 
are roughly constant. The ratio between the recombination 
coefficients slowly evolves with temperature, oc T®;® 6 ; the 
Fhi/FhoII ratio does not vary by construction if the UVB is 
spatially uniform; finally for A < 10, the gas is fully ionized 
thus nHoiii/wHii « fiHc/jiH = 0.785 is roughly constant. At 
higher densities nHeiu/nmi at the photoionization equilib- 
rium starts to decrease because He n recombines faster than 
H 11 and THeii/n-He is roughly one order of magnitude lower 
than Fm/riH- The mean rj value from our OT simulation is 
rj » 50, consistent with the expected (from eq. Illfl value, 
rj ~ 0.43 x Thi/Thch « 52. 

RT effects are crucial to detexrmine the actual values of 
rj in the IGM. This can be seen by comparing the previous 
results for the OT case with the ones from full RT calcu- 
lations, shown in the same figure; the two cases differ from 
each other both in the mean value and in shape. The RT 
mean rj value is larger than the OT one by about a factor 
of four because of the different amount of attenuation suf- 
fered by H 1 and He 11 ionizing radiation. As we have seen 
in Sec. 3, the He 11 ionizing flux is strongly suppressed with 
respect to the OT case, because even the most underdense 
voids in photoionization equilibrium with the UVB contains 
enough He 11 to remove photons with energies above 4 Ryd 
from the background radiation; the higher values we find for 
rj in the RT simulation can be interpreted in terms of the ra- 
tios between the correction factors discussed in Sec 3. The rj 
parameter shows furthermore a strong dependence on over- 
density: rj slowly grows with density, reaches a maximum 
around A ~ 7 and then rapidly declines to values which are 
still well in excess of the OT ones. In fact, at low density 
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riHem/«.Hii is roughly constant, and density dependence is 
governed by Thi/PhoU- This ratio increases with density as 
the universe is more opaque in He n than in H i , as already 
discussed. When the density is high enough, nHein/fiHii de- 
creases more rapidly than the growth of Thi/Timi, due both 
to the latter effect and to He in recombining faster than 
H n . This results in the inverted trend of the rj dependence 
on A. 

We compare our results from RT simulations with ob- 
servational data by Reimers etal. (2004) in Fig. [§] using the 
relation A = 0.8AoM,i3 (Madau, Ferrara & Rees 2001). As 
the data released are still preliminary and could suffer from 
errors associated with the analysis technique (Fechner etal. 
, 2004; private communication), we here consider the mean 
rj value for A < 7, as at higher densities the uncertainties 
on the rj values inferred from observations become very high 
and metal line contamination is a problem. Although it is 
important to point out that the dispersion in the data is 
substantial. 

The comparison in Fig. [5] clearly shows that the mean rj 
value from observations is inconsistent with the OT results, 
i.e. the expected value for the case of a spatially uniform 
UVB. The higher rj measured value can be explained when 
accounting for the radiative transfer of the metagalactic ion- 
izing radiation. Furthermore at A > 7, the RT simulation 
output predicts an 77-A anti-correlation. This suggests that 
RT effects, such as shielding and shadowing, can explain 
the systematically lower rj values observed in mildly/highly 
overdense absorption systems. According to our results, the 
lower values of rj associated with high density absorption 
systems do not arise, as previously thought, from the hard- 
ening of UVB spectrum produced by the presence of local 
hard-spectrum sources. Instead, we find that the metagalac- 
tic ionizing flux is softer in filaments than in voids. Our pre- 
ferred interpretation of the rj trend seen in Fig. [5] is that the 
rj-A anti-correlation is a direct consequence of the fact that 
in high density regions riHeii increases much faster than rim- 

3.3 Effects on the IGM Temperature 

Fluctuations in the photoionization rates could in principle 
affect the IGM temperature. In this paragraph we attempt 
to quantify such effect. The upper panel in Fig. llOl shows the 
ratio between the temperature in the RT simulations, Trt, 
and the initial temperature, Tot (estimated in the SPH sim- 
ulation in the OT limit), along a random line of sight (los); 
the lower panel shows the overdensity, A, along the same 
los. The temperature ratio Trt /Tot is significantly differ- 
ent from unity only in regions corresponding to high density 
peaks, where the temperature shows deviations of the or- 
der of 1 — 2%. The temperature in low density regions is 
basically unchanged. This is not at odd with previous stud- 
ies (Abel & Haehnelt 1999; Bolton etal. 2004) on the effect 
of radiative transfer on the temperature of the ionized gas, 
which indeed suggested that RT can be effective in boost- 
ing the IGM temperature by at least factor of two during 
the reionization process or in general when considering the 
photo-ionization of a neutral gas, with respect to the case in 
which a the gas is instantaneously photoionized and photo- 
heated by a uniform UVB. As discussed in the Introduction, 
such boost in temperature is due to the filtering of the ion- 
izing radiation. The latter is significantly harder at the edge 
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Figure 9. The parameter rj is shown as a function of the overden- 
sity A: the solid (dotted) line represents the mean rj value from 
the RT (OT) simulation and the dashed curves correspond to 
la r.m.s. deviation of the RT distribution. The long-dashed line 
shows the mean rj value from Reimers et. al. (2004) in the range 
of the observed densities (we exlcude data points with A > 7, as 
at higher densities the uncertainties on the rj values inferred from 
observations become very high and metal line contamination is a 
problem). 
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RT simulations, Trt, and the temperature estimated from the 
SPH simulation in the OT limit, Tot along a typical light of 
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of the I-front, getting harder as the I-front moves further 
from the source; this is due to the fact that the large cross- 
section, low-energy photons are selectively absorbed inside 
the HII region where they compensate for recombinations. 
The harder spectrum at the I-front yields a higher heating 
power with respect to analogous one provided by the unfil- 
tered spectrum and, due to the long characteristic cooling 
times, the low density IGM retains memory of the initial 
higher temperature phase. In the present study we focus on 
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the transfer of ionizing radiation in a ionized IGM, thus the 
above filtering process is not efficient due to the low IGM 
opacity. In order to account properly for RT effects on the 
IGM temperature, a full and consistent simulation of the 
reionization process including RT is necessary. 



4 CONCLUSIONS 

In this study we have analyzed the effects of radiative trans- 
fer on the spatial and spectral properties of the ionizing 
metagalactic radiation at a fixed redshift, z = 3.27. So far, 
this is the first attempt to study the properties of the ultravi- 
olet background radiation including full 3D RT calculations. 
Such investigation is motivated by the increasing number of 
observational evidences for a fluctuating UVB (Shull etal. 
2004; Reimers etal. 2004). These observations reveal sig- 
nificant fluctuations on scales of few comoving Mpc, which 
could hardly be originated from the inhomogeneous spatial 
distribution of ionizing sources, whose mean separation is 
« 100 h' 1 Mpc for QSOs at z w 3. 

We here neglect the UVB inhomogeneities arising 
from source clustering or non-uniform emissivity properties 
(which have been shown to be relevant only on large scales 
and at high redshift, i.e. z > 4), and focus on smaller scales. 
Our approach consists in performing RT simulations, ap- 
plying CRASH to a cosmological density field exposed to an 
ionizing background characterized by a spatially constant 
emissivity. Such radiation field, by interacting with the fila- 
ments of the cosmic web, is rapidly turned into a fluctuating 
one. 

We found that RT effects result in relevant fluctuations in 
the photo-ionization rates, whose structure is closely linked 
to the underlying density fluctuations. The UVB mean in- 
tensity, and thus its ionizing power, is progressively sup- 
pressed moving towards higher densities; such effect is found 
to be more important for radiation above 4 Ryd, due the 
high He n opacity of the universe. The ultimate physical 
nature of the decreasing trend of photoionization rates with 
density is a mixture of self-shielding and shadowing effects. 
Interestingly, we found that the STm dispersion increases 
towards higher A (from 8% at A = 0.3 to roughly 15% at 
A = 3), which is not the case for SFuc- The differences in 
H i and He n photo-ionization rates originate from the fact 
that while the H i optical depth of the IGM is dominated 
by rare, overdense regions, the He n optical depth is non- 
negligible even in voids. 

We further found that shielding in overdense regions (A > 7) 
from He n ionizing radiation, produces a decreasing trend 
of the r\ parameter with overdensity; this trend has already 
emerged from observational data (Shull etal. 2004; Reimers 
etal. 2004). According to our results, the lower values of r\ as- 
sociated with high density absorption systems do not arise, 
as previously thought, from the hardening of UVB spectrum 
produced by the presence of local hard-spectrum sources. In- 
stead, RT simulations suggest that the r]-A anti-correlation 
is a direct consequence of the fact that in high density re- 
gions riHeii increases much faster than nm- A stronger sup- 
port for the presence of RT induced UVB fluctuations comes 
from the comparison of the mean rj value from Reimers etal. 
(2004) observations and our numerical results. We find in 
fact, that the mean observed r/ value is too high to be con- 



sistent with a uniform UVB and that, accounting properly 
for the filtering of the ionizing radiation trough the translu- 
cent IGM, it is possible to reconcile the observational data 
with theoretical expectations. 

The results of the present study have several impli- 
cations. In fact, the wealth of information derived from 
the Lya forest has been so far extracted under the as- 
sumption of a uniform UVB, both from real and simulated 
data; the current results should thus be revised in order 
to account for the more realistic scenario of a fluctuating 
UVB. We give some examples. The matter power-spectrum 
inferred from the observed/simulated trasmissivity of the 
Lya forest (Croft etal. 1999) provides direct information on 
the primordial universe. UVB fluctuations could affect the 
power-spectrum shape derived in previous studies. Gnedin 
& Hamilton (2002) used RT simulations to study the effects 
of UVB fluctuations induced by quasars inhomogeneous dis- 
tribution on the matter power-spectrum (however, they did 
not look directly at the RT effects on the spatial UVB fluc- 
tuations), not finding a significant variation. However, the 
question of whether small scale fluctuations in the UVB 
could be relevant in determining the shape of the matter 
power-spectrum needs further study. 

The UVB fluctuations we have found could provide an 
alternative explanation, other than He n reionization oc- 
curring at z ?s 3, for the observed fluctuations in the r\ 
parameter as well as for the ones in the abundance ratio 
C iv /Si iv . As discussed in the Introduction, it has been 
shown that current data are not inconsistent with a scenario 
in which He n is completely reionized as early as z ~ 6 by 
thermal radiation produced in collapsing structures (Miniati 
etal. 2004); the observed fluctuating C iv /Si iv abundance 
ratio and the He n patchy opacity could be then explained 
as direct consequences of UVB fluctuations arising from ra- 
diative transfer of the metagalactic ionizing radiation trough 
the filamentary IGM. 

Although in the present study we have carried out our 
analysis at z ~ 3, we can infer that the UVB fluctuations 
become more important towards higher redshift due to the 
increasing mean opacity of the IGM. In a future work we 
plan to investigate accurately the redshift evolution of UVB 
fluctuations. 

Finally, it is worth reminding that we have neglected 
the possible contribution of local sources embedded in 
high density regions to the ionizing radiation. The rel- 
ative contribution of local sources embedded in quasar 
absorption line systems, which are commonly identified as 
high density peaks in the IGM, has been subject of recent 
debate. Schaye (2004) argued that contribution from local 
source dominates the UV metagalactic radiation in these 
absorption system, while Miralda-Escude (2004) has come 
to the opposite conclusion, deriving an upper limit to the 
importance of local sources by a consequence of the surface 
brightness conservation. The problem seems still unsettled 
and it would be desiderable to further investigate it via 
numerical RT simulations including both the UVB and the 
local sources radiation. 
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